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A CONSTANT REGRESSION CHARACTERIZATION OF A 
MARCHENKO-PASTUR LAW 

KAMIL SZPOJANKOWSKI 


Abstract. Lukacs type characterization of Marchenko-Pastur distribution in free probability is studied 
here. We prove that for free X and Y when conditional moments of order 1 and — 1 of (X + Y) —1 ' 2 X(X + Y) —1 1 
given X + Y are constant then X and Y have Marchenko-Pastur distribution. 


1. Introduction 

Since the paper [T3] free probability theory has been developed in many various directions. It turns out 
that many classical results for independent random variables such as for example Central Limit Theorem 
have their free analogues. One of the deepest known relations between classical and free probability are 
so called Bercovici - Pata bijections which give bijection between infinitely divisible distributions in free 
and classical convolution. 

In this paper we are interested in characterization problems in free probability. This seems to be 
another field which gives interesting connection between classical and free probability. Our result is a new 
example of known, but not completely well understood phenomena of analogies between characterizations 
in classical and free probability. A basic example of such analogy is the Bernstein’s theorem which 
characterizes a Gaussian distribution by independence of X + Y and X — Y for independent X and Y. 
In [7] it is proved that similar result holds for the Wigner semicircle law when independence is replaced 
by freeness assumption. 

The main result of this paper is closely related to the Lukacs theorem which provides a characterization 
of a Gamma distribution by independence of V = X + Y and U = X/{X + Y) for independent X and 
Y (see 0). It is known that the assumption of independence of U and V can be replaced by a weaker 
assumption of constancy of regressions E ([/|E)andE {U 2 \V) (see 0). In [T5J it is proved that constancy 
of regressions E (U\V) and E (t/ -1 |E) also characterizes a Gamma distribution. 

The Lukacs property was also studied in a context of free probability in T|, where Laha-Lukacs 
regression of free Meixner family is studied (see also [SH]). Theorem 3.1 from [T] contains as a special 
case a free analogue of Lukacs regressions in the case of constancy of regressions of U and U 2 given 
by V. It turns out that such conditions characterize the Marchenko-Pastur (free Poisson) distribution. 
The proof that the Marchenko-Pastur distributed X and ¥ have the property that V = X + Y and 
U = (X +Y) -1 / 2 X(X +Y) -1 / 2 are free can be found in [10]. In this paper we prove a free analogue of the 
Lukacs regressions in the case of regressions U and U -1 given by V. The proof of the main result relies 
mainly on the technique developed in our previous papers mm- 

The paper is organized as follows: in the Section 2 we briefly introduce basic notions of free probability 
and known facts which are needed to prove the main result. In Section 3 we state and prove the main 
result of the paper. 


2010 Mathematics Subject Classification. Primary: 46L54. Secondary: 62E10. 

Key words and phrases. Lukacs characterization, free Poisson distribution, free cumulants. 


1 



2 


K. SZPOJANKOWSKI 


2. Preliminaries 

In this section we give a collection of facts which we need in this paper. For more detailed introduction 
we refer to our previous papers mm- A comprehensive introduction to free probability can be found in 
Il4l or jS], 

By a non-commutative probability space we understand a pair ( A , <p) where A is a unital algebra and 
ip is a faithful, normal, tracial state. 

For self-adjoint and free X, Y with distributions /i and iq respectively, and X positive, that is the 
support of /r is a subset of (0, oo), free multiplicative convolution of p and v is defined as the distribution 
of \/X Yy/X and denoted by pMv. Due to the tracial property of <p the moments of YX, X Y and y/X Yy/X 
match. 

Let x = {-61,-82 ,...} be a partition of the set of numbers {1,..., fc|. A partition x is a crossing 
partition if there exist distinct blocks B r , B s £ \ and numbers ii,z 2 £ B r , ji,j 2 £ B s such that 
ii < j 1 < i 2 < j'2■ Otherwise x is called a non-crossing partition. The set of all non-crossing partitions 
of {1,..., k} is denoted by NC(k). 

For any k = 1,2,..., (joint) cumulants of order k of non-commutative random variables Xi,..., X„ 
are defined recursively as fc-linear maps TZk : C (aq, i = 1,..., k) —> C through equations 

¥>(Y 1; ..., Y m ) = e n 6-|s|(Yj, i £ B) 

X€NC(m) Bex 

holding for any Y,; £ {Xi,..., X„}, i = 1,..., to, and any m = 1,2,..., with \B\ denoting the number of 
elements in the block B. 

Freeness can be characterized in terms of behavior of cumulants in the following way: Consider unital 
subalgebras (Ai)i^i of an algebra A in a non-commutative probability space (A, p). Subalgebras (Ai)izi 
are freely independent iff for any n = 2,3,... and for any Xj £ A t (j) with i(J) £ I, j = 1,... ,n any 
n-cumulant 

8„(X 1 ,... ,X„) = 0 

if there exists a pair k, l £ {1 ,..., n} such that i(k) ^ i(l). 

In sequel we will use the following formula from ,5] which connects cumulants and moments for non- 
commutative random variables 

n k 

(1) y(X 1 ...X n ) = E . E. ft fc (Xi,Xi 2 ,... ,Xj fc ) n^ + 1 ...x, +1 _ l) 

k =1 l<i 2 <...<ik<n j—1 

with i i = l and ik+i = n + 1 (empty products are equal 1). 

For free random variables X and Y having distributions p. and zq respectively, the distribution of X +Y, 
denoted by /j EB iq is called free convolution of fj, and v. 

Non-commutative conditional expectation is well defined in a so called W* probability spaces, i.e. 
non-commutative probability spaces where algebra A is a von Neumann algebra. Non-commutative 
conditional expectation has many properties analogous to those of classical conditional expectation. For 
more details one can consult e.g. [12] , Here we state two of them we need in the sequel. 

Lemma 2.1. Consider a W*-probability space (A,<p). 

• IfX£A and Y £ B, where B is a von Neumann subalgebra of A, then 

(2) <p(XY) = ip(<p(X\B)Y). 

• If X, Z £ A are free then 


(3) 


y>(X|Z) = y>(X)I. 
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Now we introduce basic analytical tools used to deal with non-commutative random variables and 
their distributions. 

For a non-commutative random variable X its r-transform is defined as 

OO 

r x (z) = ^K„ +1 (X) z n . 

n —0 

In [13J it is proved that r-transform of a random variable with compact support is analytic in a neigh¬ 
bourhood of zero. From properties of cumulants it is immediate that for X and Y which are freely 
independent 


(4) 

If X has the distribution p, then often we will write r^ instead r%- 
The Cauchy transform of a probability measure v is defined as 

v(dx) 


G v (z) = f 
Jr 


Jr z-x 

Cauchy transforms and r-transforms are related by 


3f(«) > 0. 


r»( z ) + -) = z - 


(5) G v 
Finally we introduce moment generating function M% of a random variable X by 

OO 

Mx(z) = Y, ¥>(X")* n . 

It is easy to see that 

( 6 ) 

We will need the following lemma proved in |5] . 


n —1 


z = 


Lemma 2.2. Let V be compactly supported, invertible non-commutative random variable. Define C n = 
7 Z n (V'W,..., V), and C(z) = CiZ 1-1 . Then for z in neighbourhood of 0 we have 


(7) C(z) 

where r(z) is R-transform o/V. In particular, 


z + C i 
1 + zr(z) ’ 


n— 1 

(8) C 2 = l-C 1 7e 1 (Y), C n = -Y J C i n n -iW),n> 2 

2=1 


3. Main Result 


In this section we study a regressive characterization of the Marchenko-Pastur distribution which is a 
free counterpart of the characterization of the gamma distribution proved in (Ill- 

Theorem 3.1. Let (A, ip) be W*-probability space, let X, Y be non-commutative random variables in 
(. A, ip). Assume that X and Y are free, X is strictly positive, Y is positive and there exist real numbers c 
i d, such that 

(9) y>(X|X + Y) = c (X + Y) 
and 

(10) p (X -1 | X + Y) = d (X + Y) -1 . 




4 


K. SZPOJANKOWSKI 


Then X and Y have free Poisson distributions p(cA,a) and u((l — c)A,a) respectively, where A = 
and a = +++ /or some C\ > 0. 


Proof. Multiplying © and m by (X + ¥)” and applying the state to both sides of the equations, we 
obtain for n > 0 

(11) v 3 (X (X + Y) n ) = op ((X + Y) n+1 ), 

(12) ip (X _1 (X + Y)") = d ((X + Y) 71-1 ) . 

Let us define three sequences (a n ) n >_i, ( fin)n>o and ( 6 n )n>o as, 

a n = p ((X + Y) n ), fi n = (X (X + Y) n ) and S n = ip (X -1 (X + Y)") . 

We can rewrite m and m as 

(13) (3 n — CQ?77,-|_1, 

(14) S n — d cx. n — i- 

Multiplying both sides of the above equations by z n and summing over n = 0,1,... we get 

(15) B(z) = c - (A(z) - 1) , 

z 

(16) D(z) = dz (a(z) H —> 
where 

oo oo oo 

A(z) = Y, <*»*", B(z) = Y PnZ n , D{ Z ) = Y 

n =0 n —0 n=0 

Using the formula m and freeness of X and Y, for sequence a n we get 


fin — B'l 0?n 

+7?-2 ( ot n -i + a ra _2Cii + ■ ■ • + a n -i) 

+ ■ ■ • + B-n+l- 

where lZ n = lZ n (X). 

For n > 0 this gives us 

71+1 

Pn = ^ ^ 'R'k ^ ^ &ii • • • • 

k =1 ii + . •• +1/2=71+1 — k 

Using the above equations we get 

OO OO 71+1 

B ( z ) = = zU 5Z Uk “u--- <*■ 

7i=0 71=0 k= 1 ii + ---+ifc=7i+l —k 

oo oo 

= Y zk ~ llZk a h z 11 ...a ik z Zk 

k —1 Ti—k —1 ii + ...+ifc=7i+l — k 

oo oo 

= Y zk ~ 1 ' 1Zk Yl a ilZ 11 ...a ik z lk . 

k—1 771=0 ii + ...+ifc =771 


Which means that 

(17) 


B(z) = A(z)r x {zA(z)), 
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where rx(z) = !Z n +iZ n . Note that rx is the r-transform of X. 

Next we proceed similarly with the sequence fi n 

fin = C\Oi n 

+C 2 ( OL n ~i + a n -20ti + ... + a n -i) 

+ ■ ■ ■ + C n + 1 , 


where C n = 1Z 


r-i 


n — n, n I ^ 1 "■) • • • i - 


n— 1 


, dla n > 0. 


Which means that for n > 0 we have 


n+1 


fin — C k dii ■ • • 0!i k 

k=1 i\-\-.. .-\-ik=n-\-l—k 


The above equation gives us 


00 n+1 


!>(*) = 5> n a„ = l>"X;c fc E 


CX-ii . . . OLi k 


n —0 


n —0 fc=l zi + ...+ifc=n+l—/c 


= X X 


.a%uZ ' 


fc=i 


n=k —1 ii + ...+ife=n+l—fc 


= X^-‘ftX X 




fc=i 


m =0 «i + ...+£fc=ra 


This implies that 


D{z) = A(z)C(zA(z)), 


where C(z) = C n +iz n - Using lemma [2T2l we get 

(18) 


n( \ _ 4/ - t £+~) + Ci 

(z)_ W 1 + ^)+ 2 A(z))' 


Using the equations m and m , we can rewrite m and © as 

A(z)rx(zA(z)) = c - {A{z) - 1), 


^>1 + :+++) = dz +> + +)' 


zA(z)r x (zA(z)) 

Let us define an auxiliary function h(z) = zA(z)rx(zA(z)), then we can rewrite the above equations as 

(19) h(z) = c(A(z) - 1), 

zA(z) + C 1 


( 20 ) 


Mz) 


1 + h(z) 


dz ^ A(z) + ~• 


Since ft(0) = 0, then in some neighbourhood of zero we can multiply (l20l) by 1 + ft. Taking into account 
that the equation m implies C\ = da- 1 , we get 

zA 2 (z ) + A(z)Ci — zA(z)d (1 + h(z)) — C\(l + h(z)) = 0. 
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In the above equation we can replace one function A in first and second term by which follows from 
m- After simple transformations we obtain 

r2 , \ K*) = _ c(d- 1) _ 

zA(z) Ci(l — c) — zA(z)(cd — 1) 

Recall that h(z) = zA(z)rx (zA(z)). Since r is analytic in a neighbourhood of 0 and lim~_>.o zA(z) = 0 
we get 

c(d — 1) 


( 22 ) 


rx{z) = 


C i(l — c) — z(cd — 1) 

From the equation (j9]), and the assumption that X and Y are positive we get 

c = ip(x)/ip (x + Y) e (o,i). 


Similarly freeness of X and Y gives us 

d = ip (X _1 (X + Y)) =l + <p (X" 1 ) <p(Y) > 1. 

The Cauchy-Schwartz inequality implies cd > 1. 

This means that X has free Poisson distribution with parameters A = c< ^Zi and a = cAi-c) • 
Next we shall determine the distribution of Y. 

Substituting in equation (ED h from ED. we get 

A 2 (z)z(cd — 1) + A(z)(zd( 1 — c) — Ci(l — c)) + Ci(l — c) = 0. 


Since A is the moment transform of X + Y, we can use the connection between moment and Cauchy 
transforms, and after substituting 2 := l/;we obtain 

g x+y( z )~M - 1) + G x +Y(z)d(l - c) - Gx+y(z)zCi( 1 - c) + Ci(l - c) = 0. 


Now using the equation © we get the r-transform of X + Y 

rx+¥(z) = C7i(l-c)-(cd-l)z‘ 

Using GD we get 


ty{z) = 


(l-c)(d-l) 
Gi(l — c) — (cd — l)z 


Which means that Y has the free Poisson distribution with parameters A = ^ ^ and a = 


=)' 

□ 
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